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The possibility that the fundamental constants vary 
is suggested by theories unifying gravity with other in- 
teractions (see, e.g. 0, |^ <S] and review 4]). The 
analysis of quasar absorption spectra by means of the 
many- mult iplet method reveals anomalies which can be 
interpreted in terms of varying fine structure constant a 

m Q' The first indication that a might have been 
smaller at early epoch came from the analysis of magne- 
sium and iron lines 0, 0| . Later inclusion of other lines 
belonging to many different atoms and ions (Si, Cr, Ni, 
Zn, etc.) as well as many samples of data from differ- 
ent gas clouds not only confirmed the initial claim, but 
made it even stronger r\. However, there are some recent 
works in which a similar analysis indicates no variation 
of a in quasar absorption spectra d, Q ■ These works use 
the same many-multiplet method and the results of our 
calculations of the relativistic effects in atoms, but ana- 
lyze different samples of data from a different telescope. 
It is important to include as much data as possible into 
the analysis to resolve the differences, and to verify or 
discard the claim of a varying fine structure constant. 

It is natural to analyze fine structure intervals in the 
search of variation of a. Indeed, initial searches of vari- 
ation of a in quasar absorption spectra were based on 
alkali- doublet lines (alkali-doublet method) [T(l Ull 
and on the fine structure of O III However, all of 

the present evidence for varying fine structure constant 
has come from the analysis of the _EI-transition frequen- 
cies (many-multiplet method) rather than fine structure 
intervals. These frequencies are about an order of magni- 
tude more sensitive to the variation of a . However, the 
corresponding analysis is much more complicated. One 
needs to perform accurate ab initio calculations of the 
atomic structure to reveal the dependence of transition 
frequencies on the fine structure constant. We have done 
such calculations for many atoms and ions in our pre- 
vious works |l4L Il5| . In the present work we do simi- 
lar calculations for some other atoms and ions for which 
data on quasar absorption spectra are available jl^ , and 
for which corresponding calculations have not previously 
been done. 
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We use the relativistic Hartree-Fock (RHF) method 
as a starting point of our calculations. Correlations 
are included by means of configuration-interaction (CI) 
method for many valence electron atoms, or by the 
many-body perturbation theory (MBPT) and Brueckner- 
orbital method for single valence electron atoms. The de- 
pendence of the frequencies on a is revealed by varying 
a in computer codes. 

The results are presented in the form 

= t^o + qx, (I) 

where x = (a^/ap) — 1, ao is the laboratory value of the 
fine structure constant, u) and are the frequencies of 
the transition in quasar absorption spectra and in the 
laboratory, respectively, and q is the relativistic energy 
shift that comes from the calculations. Comparing the 
laboratory frequencies, loq, with those measured in the 
quasar absorption spectra, u, allows one to obtain the 
value of a billions of years ago. 

The method of calculations is described in detail in 
our early works 14, I5j. Here we only discuss the details 
specific for current calculations. 

Some atoms and ions considered in the present work 
represent open-shell (many valence electron) systems. 
Therefore, the Hartree-Fock procedure needs to be fur- 
ther specified. The natural choice is to remove all open- 
shell electrons and start the RHF calculations for the 
closed-shell core. However, this usually leads to poor 
convergence of the subsequent CI method. Better conver- 
gence can be achieved using the so called approx- 
imation in which only one valence electron is removed. 
Since we calculate not only the ground state but also ex- 
cited states of different configurations, it is convenient to 
remove the electron which changes its state in the tran- 
sition. Single-electron basis states for valence electrons 
are calculated in the V^~^ potential of the frozen-core. 

The potential corresponds to an open-shell sys- 

tem. We include the contribution of the open shells into 
the RHF potential as if they were totally filled and then 
multiply them by a weighting coefficient. Note that this 
procedure must not destroy the cancellation of the self- 
action (we would like to remind the reader that there 
is exact cancellation between direct and exchange self- 
action in the RHF equations for the closed-shell systems). 

For the CI calculations we use a B-splined single- 
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electron basis set similar to those developed by John- 
son et aL[l7l ITsL The main difference is that we 
use the open-shell RHF Hamiltonian described above to 
calculate the B-splined states. 

There are two major sources of inaccuracy in the stan- 
dard CI calculations. One is incompleteness of the ba- 
sis set for valence electrons, and another is core- valence 
correlations. We use a fitting procedure to model both 
effects. We add an extra term into a single-electron part 
of the Hamiltonian for the valence electrons: 



Here ac is the polarizability of the atomic core and a 
is a cut-off parameter that is introduced to remove the 
singularity at r = 0. We use a = (Bohr radius) and 
treat ac as a fitting parameter. The values of ac for each 
partial wave {s,p,d) are chosen to fit the experimental 
energy levels of the many-electron atom. 

The term (0) describes polarization of the atomic core 
by valence electrons. It can be considered as a semi- 
empirical approximation to the correlation interaction of 
a particular valence electron with the core. It also allows 
us to improve the convergence of the CI calculations by 
modifying the single-electron basis states. Our calcula- 
tions for rare-earth ions [2(il l2ll | have demonstrated that 
using this term allows one to obtain good accuracy of 
calculations with the minimum number of single-electron 
basis states (one in each partial wave in the cited works). 

Below we present the details and results of calculations 
for the atoms and ions considered. All transition fre- 
quencies are presented with respect to the ground state. 
Therefore we use the term "energy levels" instead. If 
a transition between excited states is needed, the corre- 
sponding relativistic energy shift q is the difference be- 
tween the level shifts {q2~4i — 92 — Qi)- 

a. Manganese (Z — 25): The ground state of Mn'^ 
is 3(i^4s ^S'a and we need to consider transitions into 
the 3(i^4s4p configuration. Earlier we also considered 
transitions to the states of the Sd^Ap configuration [T3 |. 
Since in the present work we use a different basis set, we 
have repeated calculations for this configuration in order 
to check their accuracy. 

The RHF calculations are done in the V^~^ approxi- 
mation with the configuration of external electrons. 
The 4s, 4p and higher states are calculated in the same 
V^~^ potential. We use ac = 2.05a^ for the p-wave as 
a fitting parameter in Eq. The results are presented 
in Tabled 

We have done several tests to obtain the error esti- 
mates. Fitting changes both energies and g-coefficients 
by less than 6% for all transitions, and the agreement 
with previous calculations is within 15%. We note that 
the accuracy of the fine structure splitting is significantly 
worse than this, especially for the upper levels. In some 
sense this is unsurprising: the fine structure splitting is 
much smaller than the total relativistic effect due to the 
strong cancellation between levels. Nonetheless, we may 



TABLE L Energies and relativistic energy shifts (g) for Mn+. 



State 




ujo {cm' 




q (cm ^) 






theory 


experiment 








no fittinj 


J fitted 




this work [14] 




'P° 


36356 


38424 


38366 


869 (150) 986 


3d^4p 


7 r)o 
^3 


36519 


38585 


38543 


1030 (150) 1148 


3d^4p 


7 no 

-M, 


36755 


38814 


38807 


1276 (150) 1420 


3d'^4s4p 




84092 


83363 


83256 


-3033 (450) 


3d'^4s4p 


7 no 

-fs 


84319 


83559 


83376 


-2825 (450) 


3d'^4s4p 


7 no 
^4 


84619 


83818 


83529 


-2556 (450) 



try to fit the splitting using some parameter, and see how 
q changes. In fact we have done this by varying a itself. 
The splitting for the higher levels becomes very close to 
experiment at approximately x — —0.5, and at this point 
q is within 1% of the value at a; = 0. This is a similar 
error estimation procedure as that used in jl5l | for Ni II, 
and it shows that q is insensitive to the fine structure 
splitting. In another test, we have fit the fine structure 
splitting by introducing screening parameters before the 
Coulomb integrals in the CI code. This models the effect 
of screening of the valence electrons by the core electrons. 
Again, q changes by less than 15%. 

More worrying is the possibility of level misidentifica- 
tion. This can occur due to misidentification in the ta- 
bles, or alternatively in the CI results. The Mn+ spectra 
contains two 3d'*4s4p ''P° multiplets in close proximity. 
It is possible (although unlikely) that they are swapped 
in the computer code. The q values for the other multi- 
plet are approximately 15% smaller than for the required 
transitions. Therefore, we use 15% as a conservative es- 
timate of the accuracy of q. 

Note that the relativistic shift is positive for the s — p 
singe-electron transitions and negative for the d — p tran- 
sitions. Having transitions with different signs of q- 
coefficients in the same atom (ion) helps to fight sys- 
tematic errors in the search for variation of a 14] . 

b. Titanium( Z = 22): We perform calculations for 
both Ti+ and Ti^+ starting from the same RHF approxi- 
mation, and using the same single-electron basis set. The 
ground state of Ti+ is 3d^4s ^^3/2 and we need to con- 
sider transitions into states of the 3d^4p and 3(i4s4p con- 
figurations. The ground state of Ti^+ is 3(i^ and we 
need to consider transitions into the states of the idAp 
configuration. Therefore it is convenient to do the RHF 
calculations for the Ti-^+ ion with the 3d^ open-shell con- 
figuration. The 4s, 4p and other basis states for the CI 
method are calculated in the frozen-core field of Ti^+ . 

The fitting parameters chosen are ac — 0.38a'g for s- 
electrons and ac = 0.065a^ for d-electrons. The results 
are presented in Table ^1 As in the case of Mn"*" , there 
are negative and positive relativistic shifts. The effects 
of fitting and change of basis set does not exceed 10%. 
The values of the g-coefficients for titanium are consistent 
with calculations for other atoms and with semi-empirical 
estimations using the formulas presented in |l4j . In par- 



ticular, the values of the negative ^-coefficients for the 
d — p transitions are v ery close to the values for simi- 
lar transitions in Cr II The positive coefficients for 
Ti+ are very close to those for Mn+ after rescaling by 
according to the semi-empirical formula The 
q-coefficients for the 3cPAp '^F° states closely match in- 
dependent calculations of [23| . 



TABLE II: Energies and relativistic energy shifts {q) for Ti"* 
and Ti2+ . 



State 


LOO (cm" 




q (cm ^) 




theory 


experiment 






no fitting fitted 


[22] 





3d^4p 

3d4s4p 

3d4s4p 

3d4p 



'-^5/2 

4 T70 
4 r.o 

4 r)0 

^1/2 
4 p)0 

^3/2 
4 rjo 

^1/2 
4 po 

-^^3/2 



28097 
29401 
29521 
31143 
31227 
50889 
51341 

80558 



29759 
30691 
30813 
32416 
32510 
52185 



29544 
30837 
30959 
32532 
32603 
52339 
52330 



Ti 



2+ 



396 (50) 
541 (50) 
673 (50) 
677 (50) 
791 (50) 
-1564 (150) 
-1783 (300) 



77000 -1644 (150) 



Figure ^ shows the dependence of the highly excited 
3d4s4p levels of Ti+ on . The data points are gen- 
erated from a highly saturated CI, without any fitting 
parameters. The dependence is very close to linear, 
and hence there is no level pseudo-crossing in the code. 
This does not completely exclude the possibility that the 
^F°^2 ''£'3/2 levels are strongly mixed (this could be 
checked and accounted for if the experimental g-factors 
were known but they are not available). In the case 
of strong mixing, q lies somewhere between that of the 
two states, thus we have increased the '^F^i2 error to in- 
clude this possibility. 

c. Sodium (Z —11): In contrast to the ions consid- 
ered above, sodium has only one external electron above 
closed shells. Its ground state is ls^2s^2p^3s ^5i/2- Very 
accurate calculations are possible for such systems by in- 
cluding certai n ty pes of correlation diagrams to all orders 
(see, e.g. 0,l23)- However, since both relativistic and 
correlation effects for Na are small we use a simplified 
approach. We calculate the correlation potential S (the 
average value of this operator is the correlation correc- 
tion to the energy of the external electron) in the second 
order of MBPT, and use it to modify the RHF equations 
for the valence electron and to calculate the so-called 
Brueckner-orbitals. Note that due to iterations of S cer- 
tain types of correlation diagrams are still included in all 
orders in this procedure. The final accuracy of the energy 
is better than 1%, and for the fine structure accuracy is 
2-6% (see Table Unl). We believe that the accuracy for 
the relativistic shifts q is on the same level. 

d. Carbon (Z ~ Q): Relativistic effects for carbon 
and its ions are small and calculations can be done with- 
out fitting parameters. The ground state of neutral car- 
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FIG. 1: Dependence of some odd 3d4s4p levels of Ti+ on 
X — (a/ao)^ — 1. Solid lines correspond to J = 3/2 and 
dashed lines to J = 5/2. Data points are shown from sim- 
ulations, along with a linear fit that shows that there is no 
interaction between the transitions. The assignment of the 
multiplets, from the bottom up: "7?°, ^F°, and the 
*_D°, '^F° multiplets appear in the wrong order compared to 
experiment. The *i3i/2 hue is not shown since there are no 
nearby lines with the same angular momentum, and hence 
no risk of interaction. The ^1)7/2 and '^Fy2 9/2 lines are also 
omitted. 



TABLE III: Energies and relativistic energy shifts {q) for Na. 





State 


theory 


loq (cm ^) 

experiment ^2] 


q (cm ^) 


3p 


2 po 

'^1/2 


16858 


16956 


45 (4) 


3p 


2 po 

-^3/2 


16876 


16973 


63 (4) 


4p 


2 po 
^1/2 


30124 


30267 


53 (4) 


4p 


2 po 
^■i/2 


30130 


30273 


59 (4) 



bon is ls^2s^2p^ ^Pq. Our RHF calculations for this 
atom include all electrons, however, since we need to con- 
sider configurations with excitations from both 2s and 2p 
states, we treat both as valence states in CI. 

For neutral carbon we have performed the calcula- 
tions for the ground state configuration as well as for 
excited configurations 2s^2p3s, 2s2p^, 2s^2p4s,2s^2p3d, 
2s^2p4:d, 2s^2p5d and 2s^2p6d. However, we present in 
Table IIVI only results for the 2s2p^ configuration. The 
relativistic energy shift for all other configurations is 
small (|g| < 100 cm~^). This is smaller than the un- 
certainty in the g-coefficients for heavier ions. Since the 
analysis of quasar spectra is based on comparison of the 
relativistic effects in light and heavy atoms (ions), small 



4 



TABLE IV: Ener; 
carbon atom and 



pes and relativistic energy shifts (q) for the 
its ions (cm^^) 



State 


LJo (cm ^) q (cm ^) 




theory experiment [221 



2s2p^ 
2s2p^ 
2s2p^ 
2s2p3 
2s2p=* 

2s^2p 
2s2p2 
2s2p2 
2s2p2 

2s2p 

2p 
2p 



3po 

2 po 
-f^3/2 

Ji35/2 

D3/2 

Ipo 



2 po 

M/2 
2 po 

-^3/2 



c 



66722 
66712 
66716 
75978 
100170 

74 
76506 
76503 
97993 

103955 

64557 
64686 



C+ 

^2+ 
^3+ 



64087 


151 


(60) 


64090 


141 


(60) 


64091 


145 


(60) 


75254 


111 


(60) 


1 nt:;7QQ 


1 

iOU 


(bU) 


63 


63 


(1) 


74930 


179 


(20) 


74933 


176 


(20) 


96494 


161 


(30) 


102352 


162 


(20) 


64484 


104 


(20) 


64592 


232 


(20) 



tential of the closed-shell (helium) core. As can be seen 
from Table Hvl accuracy for the energies is within 10%. 
We estimate the accuracy of g-coefficients at around 10%. 



TABLE V: Energies and relativistic energy shifts (g) for oxy- 
gen ions. 



State 


theory 


tJo (cm ^) 

experiment [2^ 


CI ( C'vr\ ^ ^ 








0+ 






2s2p^ 




122620 




119873 


346 (50) 


2s2p'' 




122763 




120000 


489 (50) 


2s2p'* 


'Pl/2 


122848 




120083 


574 (50) 






02+ 




2s2p^ 




121299 




120058 


723 (50) 


2s2p^ 


3po 


143483 




142382 


726 (50) 








03+ 




2s2p^ 


^P>3I2 


129206 




126950 


840 (50) 






05+ 






ls^2p 


2 r>o 
^1/2 


96501 




96375 


309 (50) 


ls^2p 


2 po 

-^3/2 


97091 




96908 


913 (50) 



relativistic energy shifts in light atoms can be neglected. 
The q-coefficients for the 2s2p^ configuration are larger 
because this configuration corresponds to the 2s — 2p 
transition from the ground state. These are the lowest 
valence single-electron states with the largest relativistic 
effects. Other excited configurations correspond to the 
2p — ns or 2p — nd (n > 3) transitions. Ifowever, rela- 
tivistic energy shifts for higher states are smaller jTj| . 
The calculations for C^"'' and 0"^+ are done in the po- 



e. Oxygen (Z = 8): Relativistic effects for oxygen 
ions are comparatively large, and become larger with in- 
creasing electric charge. This is in agreement with semi- 
empirical formulae presented in 14]. For transitions in 
neutral oxygen, however, \q\ < 20 cm~^; these results are 
not presented here. 
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